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Abstract. The first paper attempting direct calculation of the Yang-Mills vacuum wave 
functional was published by Greensite in 1979. I review some recent results of the determination 
of the vacuum wave functional in Monte Carlo simulations of SU(2) lattice gauge theory. 


1. Introduction 

The DISCRETE 2014 Symposium commemorated a number of important events that had 
shaped modern physics. Maxwell presented his theory of electromagnetism to the Royal Society 
150 years ago; non-abelian gauge theories were proposed by Yang and Mills 60 years ago; the 
Brout-Englert-Higgs mechanism, the quark model, and Bell inequalities are 50 years old; CP 
violation was experimentally discovered 50 years ago as well. In my talk I would like to draw 
attention to another anniversary: In February 1979, the first paper trying to calculate the 
ground-state wave functional of Yang-Mills (YM) theory was submitted to Nuclear Physics 
by Jeff Greensite [lj. Thirty-five years have passed since then, but the problem still defies 
satisfactory solution. 

Formulation of the problem The vacuum wave functional (VWF) 'I'o of quantum 
chromodynamics in the Schrodinger representation depends on quark fields of six flavours with 
three colours, each represented by a Dirac four-component bispinor, and on eight four-vector 
gluon fields - this is altogether 104 fields at each point in space (not taking constraints from 
gauge invariance into account). This is a formidable object from both mathematical and practical 
point of view. To simplify the problem, one can reduce the number of colours from three to 
two, omit quarks, discretize space ( i.e . formulate the theory on a lattice), and eventually go to 
lower-dimensional spacetime. One can hope |2j that the resulting model captures at least gross 
features of the full theory, in particular information on the mechanism of colour confinement. 

Omitting quarks, the SU(2) YM Schrodinger equation in (d + 1) dimensions in temporal 


* Based on talks presented at the Fourth Symposium on Prospects in the Physics of Discrete Symmetries, 
DISCRETE 2014, King’s College, London, UK, December 2-6, 2014, and (under the title Measurement of the 
Yang-Mills vacuum wave functional in lattice simulations) at the fth Winter Workshop on Non-Perturbative 
Quantum Field Theory , INLN, Sophia Antipolis, France, February 2-5, 2015. Abridged versions will be published 
in proceedings of these conferences. 


gauge looks very simple: 
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Physical states are simultaneously required to satisfy Gaufi’ law: 

(5 ac d k + ge abc A\) ^S>[A\ = 0. 
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A few well-known facts If we set the gauge coupling g to 0, the Schrodinger equation 0 
reduces to that of (three copies of) electrodynamics and its ground-state solution is known to 

be [3]: 
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Due to ([2]), the ground state must be gauge-invariant. The simplest form consistent with ([ 3 ]) is 
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where A is some adjoint-representation kernel that reduces, in the limit g —> 0, to (—A) x / 2 . 

At long-distance scales, one expects the VWF to be the state of magnetic disorder B0U5]: 

*o[^]~Afexp [~\hj t^x F?j(x) F$(x) }. (5) 

This is also called the dimensional-reduction (DR) form. With such a VWF, the computation 
of a spacelike loop in (d + 1) dimensions reduces to the calculation of a Wilson loop in YM 
theory in d (euclidean) dimensions. If the vacuum were of the DR form for YM theories in both 
(3 + 1) and (2 + 1) dimensions, then these would be confining, since the theory in 2 euclidean 
dimensions exhibits the area law. However, this cannot be the whole truth, for various reasons. 
Via DR, one gets e.g. the area law for colour charges from all representations r of the gauge 
group, and string tensions a r of the corresponding potentials are proportional to eigenvalues 
of the r-th quadratic Casimir operator. Approximate Casimir scaling is observed at short and 
intermediate distances, but at large distances, due to colour screening , string tensions depend 
on the JV-ality of the representation (see figure [I] for illustration). 


Approaches to the problem A number of different strategies were adopted in attempts to 
determine the YM vacuum wave functional (see |6j for a more extensive set of references): 

(i) Strong-coupling expansion of the VWF mm- 

(ii) Weak-coupling expansion of the VWF mm- 

(iii) Formulation of the theory in cleverly selected variables and expansion of the VWF in terms 
of these new variables j!3l4l9j . 

(iv) Variational Ansatz for the VWF in a certain gauge; determination of its parameters by 
minimizing the expectation value of the YM hamiltonian (see e.g. mm)- 

(v) Guess of an approximate form of the VWF and tests of its consequences (see {22j, [23U25] ). 












Static quark-antiquark potential, SU(2) gauge theory 



Figure 1. Sketch of static potentials in 
SU(2): At short distances the potentials are 
Coulomb-like. At large distances, colour 
charges of higher-representation sources can 
be screened by gluons. Potentials of integer 
representations become asymptotically flat, 
while asymptotic string tensions for all 
representations with half-integer j are the 
same as for j = 1/2. Casimir scaling of string 
tensions is expected only at intermediate 
distances, at least in some models of 
confinement. 


2. Flatland: A romance of two dimension^] 

I will start with the problem in the simplest setting: with SU(2) Yang-Mills theory in “ Flatland ”, 
i.e. in three spacetime dimensions. That model deserves attention not only because it is more 
tractable than the realistic case, but also due to its relevance to the high-temperature phase of 
chromodynamics in four dimensions^ 


Hints on the form of the VWF In the lattice formulation, a systematic strong coupling 
expansion of the YM VWF is of the form ’Fo = A/" exp(—i? [[/]), where the function R in the 
exponent is an expansion in terms of closed loops, products of link matrices U along closed 
contours on the lattice [7J- Guo, Chen and Li [8j computed the first few terms of this expansion 
and showed that for slowly varying fields they organize themselves into the following serief0 

R[U] oc W)Tr [B 2 ] - ^Tr [B(-P 2 )B] + ... (6) 


Here go and g\ are functions of the lattice spacing a and the coupling constant g, B = F 12 is the 
colour magnetic field strength, and V 2 = T> j. • T>\. is the adjoint covariant laplacian, where T>y [A] 
denotes the covariant derivative in the adjoint representation. The first term of ([b]) corresponds 
to the dimensional-reduction vacuum wave functional ([5]). One can imagine that an expansion 
of the form ([b]) might come from the wave functional 
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with the kernel introduced in Q, being a functional of the adjoint covariant laplacian. 

A similar hint may be deduced also from the approach of Karabali et al. |13| . They combine 
two components of the gauge potential into complex-valued fields {A, A} = ^ (Ai ± iAf) and 
introduce new variables, a matrix-valued field M e SL (N,C), related to A, A via 


A = — (cbM) M 1 , A = 1 ^cblVF^ , where {z, z} = x\ ± 1 x 2 - (8) 


1 This section’s heading paraphrases the title of a novella [26] by English schoolmaster and theologian Edwin 
Abbott Abbott, writing pseudonymously as A Square , which describes a two-dimensional world inhabited by 
geometric figures. The book was published in 1884 by a publishing house at that time located in Essex Street, 
less than half a kilometre from King’s College. 

2 Another good justification - not scientific, but utterly human - for investigating the model in lower dimensions 
was given by Feynman in his lecture at the 1981 EPS-HEP Conference in Lisbon m-- “Of course, understanding 
something in 2+1 dimensions does not imply that you understand it in 3+1 dimensions, but, after all I have 
lots of time, I have tenure, so I can do whatever I want. ” I would avoid writing such a justification into a grant 
proposal, but find it nevertheless very reasonable. 

3 The trace in this symbolic expression includes sums over colour indices and lattice sites. 





















Figure 2. Ao vs. (3 from simulations of 
3D euclidean SU(2) lattice gauge theory at 
various couplings /3 and lattice volumes L 3 . 
The best fit to data is Ao oc /3 -1 - 44 , which 
differs from the expected f3~ 2 dependence and 
indicates that Ao diverges in the continuum 
limit. 


M transforms covariantly, M —> J2M, under a gauge transformation fl, and is used to define 
gauge-invariant variables: 

H = M t M, J a = Tr (r a (<9,H)H“ 1 ) , (9) 

through which one expresses the hamiltonian, inner products of physical states, and the VWF. 

Karabali et al. further show that the part of the VWF bilinear in variables J a , when expressed 
through original colour magnetic fields, takes on the form: 
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( 10 ) 


This expression, however, is not gauge invariant, but one can assume that higher-order terms 
in J a might turn the ordinary laplacian in (10) into the adjoint covariant laplacian, and thus 


convert (10) again into the form (7) with R = (^/—V 2 + m 2 + m 2 ^j 


-l 


The proposal of Samuel Almost 20 years ago, Samuel [22] put forward a simple vacuum 
wave functional of the form Q, interpolating between weak-coupling ([3]) and DR ([5]) limits: 


T r o[A] = A/" exp — \ 
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He estimated with its use the 0 ++ glueball mass to be about 1.5 GeV. However, this particular 
form may be flawed if taken at face value: there are hints that the adjoint covariant laplacian 
needs to be regularized. We have computed its eigenvalues in numerical simulations of the three- 
dimensional euclidean SU(2) YM theory on a lattice. (— T> 2 ) has a positive definite spectrum, 
finite with a lattice regulator. If its lowest eigenvalue Ao were finite in the continuum limit, it 
should scale, as a function of /3 = 4/g 2 , as /3~ 2 for large (3. Our data (see figure^ indicate that 
liing^oo /3 2 Aq(/3) —> oo, i.e. Aq diverges for typical configurations in the continuum limit. 


The GO proposal We proposed [23] another simple form with a seemingly small, but crucial 
difference from that of Samuel: 
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As a remedy to the problem mentioned above, it is suggested here to subtract from the covariant 
laplacian its lowest eigenvalue. The proposed VWF contains a single free parameter, a mass m, 
that should vanish in the free-field limit (for g —> 0). The expression is assumed to be regularized 
on a lattice, and we use the simplest discretized form of the adjoint covariant laplacian: 
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where 


Uf(x) = o Tr\a a U k (x)a b ul(x) 


(14) 


and Uk{x) are the link matrices in the fundamental representation. 

We have provided a number of (semi)analytic arguments in favour of the proposed VWF ( |12[ ): 

(i) T'o reproduces the VWF of electrodynamics ([3]) in the free-held limit (for g —> 0). 

(ii) The proposed form is a good approximation to the true vacuum also for strong fields 
constant in space and varying only in time [23] . Indeed, if we put such a physical system 
into a finite volume V, its lagrangian and hamiltonian are: 
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where & = \ A± x A 2 \. It is natural to look for the ground-state solution of the Schrodinger 
equation RAf = EAf in the form of 1/V expansion: 

T'o = exp[—Vi? 0 — 4?i — V~ 1 R 2 — •••]• (16) 


The leading term has to satisfy 
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where £ = y A\ ■ A\ + A 2 ■ A 2 . The colour vectors A\ and Ao dehne a plane in colour space, 
we can choose e.g. A\ = (Ai,0,0), A 2 = (A 2 cos 0, .4.2 sin0, 0), then 6 = |^4i^42 sin0|. We 
assume ^. 1,^2 to be of the same order O(A). The contribution of Ro to (16) will be 
non-negligible if VRq ~ 0(1), i.e. Rq ~ 0( 1/4"). 

On the other hand, for strong enough fields, ~ ge abc A b k and 
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where 
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Rayleigh-Ritz Calculation, (3=6, L=16 
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Figure 3. Vacuum energy density vs. mass 
parameter m from numerical simulations on 
a lattice with L = 16 at coupling /3 = 6 . 
The minimum is away from zero, at roughly 
m = 0.3. This gives a string tension which is 
a little low for /3 = 6 , but the disagreement 
should not be taken too seriously, because 
the estimate for vacuum energy ( 22 ) is only 
approximate. 


Figure 4. Mass gaps, extracted from 
equal-time correlators of colour magnetic 
fields [23], measured in recursion lattices at 
various lattice couplings, compared to the 
0 ++ glueball masses in (2 + 1 ) dimensions 
obtained via standard lattice Monte Carlo 
methods (from [28], denoted “MC”). They 
agree within a few (< 6 ) per cent. 


The leading term in our VWF (|12[) then is 


T r o = AAexp 



6 2 

\/C3 - Cl + ( rn/g ) 2 


( 21 ) 


where C and £3 are respectively the smallest and largest eigenvalues of the above matrix Z. 
For strong constant fields (\gA\ m. Ao), both ( m/g ) 2 and Ci ~ Ro ~ 0(1/1/) are negligible 
w.r.t. C 3 = £ 2 , and the expression (21) agrees with (16) and (18). 


(iii) If we split the colour magnetic field strength B(x) into “fast” and “slow” components, the 
part of the VWF that depends on B s i ow reduces to the magnetic-disorder (DR) form ([ 5 ]). 
The fundamental string tension is then easily computed as op = 3mg 2 /16 = 3m/4/3. Non¬ 
zero value of the parameter m then implies non-zero crp, i.e. confinement of fundamental- 
representation colour charges. 


(iv) One can take the mass m in the wave functional as a free variational parameter and compute 
(approximately) the expectation value of the YM hamiltonian (see section VI of |23| for 
details). The result, expressed as a sum over eigenvalues of the adjoint covariant laplacian 
in the thermalized gauge field A, is: 


(K) 



A 0 + m 2 + 


1 Ao — m 2 ^ 
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( 22 ) 


In the abelian free-held case, the optimal value of m equals Ao and Ao goes to 0 in the 
continuum limit, so the theory is non-confining. In the non-abelian case one can evaluate 
the energy density {T-L)/L 2 numerically in Monte Carlo simulations, and finds that a non¬ 
zero (finite) value of m is energetically preferred. A typical result is displayed in figure [3} 
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Figure 5. The Coulomb-gauge ghost pro¬ 
pagator at f3 = 9 on 32 2 lattice. 


Figure 6. The colour-Coulomb potential at 
f3 = 6 on 24 2 lattice, computed from con¬ 
figurations with a cut |V(0)| < 100. 


Some numerical evidence All above arguments are encouraging, but quantitative tests are 
needed to gain more confidence in the proposed form of the VWF. Such tests are provided by 
numerical lattice simulations. For readers not familiar with the method, basics is summarized 
Appendix A 


m 


We have compared a set of physical quantities computed in two ensembles of lattice gauge- 
field configurations: 


I. 


Monte Carlo lattices: Ensemble of two-dimensional slices of configurations generated by 
Monte Carlo simulations of t hree-dimensional euclidean SU(2) lattice gauge theory with 
standard Wilson action (A.4) at a coupling f3 = 4/g 2 ; from each configuration, only one 


II. 


(random) slice at fixed euclidean time was taken. These configurations are distributed with 
the weight proportional to the square of the true VWF of the theory, | , F true [[y]| 2 . 

“Recursion” lattices: Ensemble of independent two-dimensional lattice configurations 
generated with the probability distribution given by the (square of the) VWF ( |12| ) , with m 
and g 2 fixed to get the correct value of the fundamental string tension. These configurations 
can be generated efficiently by the recursion method proposed (and described in detail) 
. Essential points of the method are sketched in Appendix B| 


m 


We computed in both ensembles the mass gap [23], the Coulomb-gauge ghost propagator, and 
the colour Coulomb potential [211 ]. The latter two quantities are of particular interest because 
of their role in the so-called Gribov-Zwanziger mechanism of confinement [29l 130 ], Example 
results are shown in figures [4j [5j and [6j The agreement between quantities measured in MC and 
recursion ensembles is very reasonable j^] 

One can also measure ( i.e . compute in numerical simulations) amplitudes of various sets of 
test configurations in the true YM vacuum, and compare them with predictions based on the 
proposed VWF ©. The method for computing relative weights of configurations was proposed 
by Greensite and Iwasaki long ago [31] and will be described in section [3j The technique was 
applied in (2 + 1) dimensions to sets of abelian plane wave configurations of varying amplitude 
and wavelength, and sets of non-abelian constant configurations [6j. The obtained data agree 


4 The agreement in the case of the colour Coulomb potential is less satisfactory and we attribute it to the 
existence of exceptional configurations that are extremely difficult to fix to the Coulomb gauge. Differences in the 
measured colour Coulomb potentials in two ensembles are small, if one ensures approximately equal population 
of exceptional configurations in both of them. (See Section IV of [24] for details on this point.) 


























with expectations based on (12), but similar agreement was obtained also for some other forms 
of wave functionals which simplify to a DR form at large scales (see El)- 

Summarizing this section, we found analytic and numerical evidence that in D = 2 + 1 
dimensions our Ansatz for the VWF seems a fairly good approximation to the true ground state 
of the theory. However, we do not live in Flatland , our world is Spaceland, so I will switch for 
the rest of this contribution to pure Yang-Mills theory in four spacetime dimensions. 


3. Version in 3D or SpacelancQ 

The extension of the proposed vacuum wave functional ( |12|) to (3 + 1) dimensions is straight¬ 
forward: one replaces the product B a (x) B b (y ) by ^F-j(x) F b j(y), and two-dimensional integrals 
by three-dimensional ones: 
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Of course, the adjoint covariant laplacian now involves ^-summation over three space directions: 
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(24) 


Positive news, but. . . The wave functional (23) is also in (3 + 1) dimensions exact in the 
free-field limit. It also approximately solves the YM Schrodinger equation in the zero-mode, 
strong-field limit. It is now convenient to introduce: 
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Then the leading term Rq of the 1 /V expansion (16) of the exponent of T'o is again given by 
since 
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(28) 


for strong fields from the “abelian valley” where the components Ai,A 2 ,A% are nearly aligned, 
or antialigned, in colour space. 


In this same limit, the proposed VWF (23) reduces to T'o = AAexp(— Q) with 

/ § ab 


Q « \gV{Ai x Aj) a 
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(29) 


The first term dominates and gives VRq ~ 0(1) as above, the other two are of order O (1/V) 
(see Appendix A of [23]). 

However, the transition to three space dimensions causes complications associated with the 
Bianchi constraint that colour magnetic fields F“- have to satisfy. Because of that constraint, 
direct numerical generation of gauge-held configurations with the probability distribution given 
by the square of the vacuum wave functional (23) is much more challenging. I have not yet been 
able to find a way of generating such configurations, similar to the recursion method used in two 
space dimensions (see section [2] and Appendix B). A different approach to testing our proposal 
is required. 

5 Not surprisingly, one can find also a book with Spaceland in the title (32]. It is a science-fiction novel written 
by American mathematician and computer scientist Rudy Rucker as a tribute to Edwin Abbott’s Flatland. 

























Morse potential, V(x)=V 0 (e' 2x -2e' x ), V 0 =2 



x 


Figure 7. Result of a numerical determi¬ 
nation of \ipo(x)\ 2 for the Morse potential 
V(x) = V 0 (e~ 2x - 2e~ x ) with V 0 = 2 [3S]. 
The red solid line shows the exact result, the 
green points are the outcome of a Monte Carlo 
computation based on (32), and dotted lines 
indicate the ±10% band around the exact so¬ 
lution. 


Measurement of the ground-state wave function in quantum mechanics Let us start 
with a simple question: How can one numerically compute the wave function of a ground-state 
of a one-dimensional system with hamiltonian H in quantum mechanics? The most direct 
approach is to solve the Schrodinger equation numerically. A less accurate way, but useful for 
our purposes, is to start from a simple formula: 


[^(x)] 2 = 1™ e E ° T G(x, 


—It; x, 0 ) = lim 

T—>• OO 


G(x, —ir; x, 0) 
fdtG(t,-i T ;t, 0 )’ 


(30) 


where G(x 2 , —ir;x 1 ,0) is the Green’s function (propagator) of a free particle moving from the 
point x\ at t = 0 to the point X2 at euclidean time (— ir). Expressing the ratio in (30) through 
path integrals in the usual way (see e.g. m) 


G(xn = x, — ir; xq = x, 0) / dx\... dxN-i exp [— J 0 r Hdr'\ 

f dx G(xn = x,—ir;xo = x,0) f dx\ ... dx^-idxN exp [— / Q r Hdr'} 

one finally gets 

\M*)\ 2 = \Jmt)} 5 E(O) - x] e- s N , (32) 

where S is the euclidean action corresponding to the hamiltonian H . This equation expresses 
the wave function squared as an average of a (5-function over all paths, a procedure that might 
appear totally inappropriate for numerical computation. However, with a simple trick it can be 
implemented efficiently [33], and a sample result is displayed in figure [7j 


The relative-weight method The method of Greensite and Iwasaki m is based on a 
generalization of ( |32| ) to quantum field theory. The squared VWF of the pure YM theory 
is given by the path integral} 


*0 P'] = \j[DU] n«x,0) - E/(X)] exp (-£[[/]), 


X,Z 


(33) 


The relative-weight method EH enables one to compute ratios T r Q[17^]/T , g[L r ^ m )] for 
configurations belonging to a finite set U = j (x), j = 1,2 ,..., M j using a simple procedure: 
One performs Monte Carlo simulations with the usual update algorithm (e.g. heat-bath) for all 

6 One could insert into the path integral a factor imposing fixing to lattice temporal gauge, which sets all timelike 
links to 1 except on one time slice at t ± 0. However, this gauge fixing is in fact not necessary. 
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Non-abelian constant conf’ns, p=2.2, 20 4 lattice 
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Non-abelian constant conf'ns, p=2.4, 24 4 lattice 
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Figure 8. [~\og(N n /N tot )] vs. R n = /iftn 
for non-abelian constant configurations with 
k = 0.14; /r at (5 = 2.2 on 20 4 lattice was 
4.06(4). 


Figure 9. [-log(N n /N tot )] vs. x n = un for 
NAC configurations with /c = 0.14 at /3 = 2.4 
on 24 4 lattice; the slope g, from a linear fit 
comes out 6.04(2). 


spacelike links at t 0 and for timelike links. Once in a while one updates the spacelike links 
at t = 0 all at once selecting one configuration from the set U at random, and accepts/rejects it 
via the Metropolis prescription. Then, for a large number of updates N to t, 


*2[C/M] W t ot~>oo N m ’ 


(34) 


where N m ( N n ) denotes the number of times the m-th (n-th) configuration is accepted. To ensure 
a non-negligible acceptance rate of Metropolis updates for all configurations in the set U , they 
must lie close in configuration space. This limits somewhat the applicability of the method. 

If the VWF is assumed to be of the form T , o[C r ] = ACexp(— R[U]), then the measured values 
of [—log(A r n /A ? tot)] should fall on a straight line with unit slope as function of R n = R[U^]. 
An example is shown in figure [8] for a set of non-abelian constant configurations (to be specified 
below). 


Direct measurement of the vacuum wave functional We used the relative-weight 
method to calculate the squared WVF for the following two classes of lattice gauge-field 
configurations: 

I. Non-abelian constant (NAC) configurations: 

C/nac = | ul n \x) = \Jl- (u( n )) 2 1 + m (n) er fe j , (35) 

where 

ne {1,2,..., 10}. (36) 

The constant k, regulating amplitudes of NAC configurations, is selected so that the ratio of 
the smallest to the largest weight within the set is not too small, at most O (10~ 4 4- 10~ 3 ). 
The expected dependence of relative weights on n and n is linear: 


= ( 
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(37) 
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Abelian plane waves, p=2.5, 30 4 lattice, n=(0,1,5) 


* 20 lattice / 



3 

-1-1-1-1-1- 



B / 




data 


10 

- AF / 



2.8 

0.975+2.07 x 






75 2.6 



=1 

a. . ■. / 



Z 

ZT 2.4 







CT> 







° 2.2 







2 







1.8 

, , , , , 


0.1 

-1-1-1-1-1- 




0.5 0.6 0.7 0.8 0.9 


3 0.5 1 1.5 2 2.5 

3 



Xj=0.5(a+Yj), a=0.73, y=0.11 


P 





Figure 10. The quantity /i, extracted from 
computed weights of NAC configurations, vs. 
the coupling (3. 


Figure 11. An example plot of 

[- log(A4 J) /iVtot)] vs. \{a n +7raj) for APW 
configurations. 


The data are consistent with this expectation; the constant ^ at a given coupling fi is 
obtained as the slope of a linear fit of [— log(N n /Ntot)] vs. nn, see figure [9j This constant 
coincides with the parameter /x that appears in the dimensional-reduction Ansatz for the 
VWF 0. 

Predictions for NAC configurations resulting from the DR and our proposed vacuum wave 
functional are identical (up to parameter-naming conventions). Moreover, the dependence 
of n on j3 was already computed in the pioneering work of Greensite and Iwasaki m 
more than 25 years ago, albeit at rather small-size lattices, 6 4 and 8 4 . New results from 
simulations on 20 4 and 24 4 lattices confirm their findings, cf. figure 2 of [31] with our 
figure 10 The strong-coupling prediction /x(/3) = 0 is confirmed at small /3, and in the 


weak-coupling region /x(/3) behaves as a physical quantity with the dimension of inverse 


mass: 


where 


rtP)f(P) = P Ph ys ~ 0.0269(3), 


m = 



exp 



(38) 


(39) 


II. Abelian plane-wave (APW) configurations: 


Uapw = S Ui ^ (x) = \l 1 - 


( w^\x) 


1 +iwil\x)cr 3 , U 2 J> (x) = U^’(x) = l}, (40) 


r(i) 


(J)/ 


where 


t»g> = 


'«n +7 nj 


LA 


COS 


27r 

T 


-n ■ x 


n = (ni,n 2 ,n 3 ), 


j e {1,2,..., 10}. (41) 


Amplitudes in a particular set of plane waves with wave number n are parametrized by 
a pair (a n ,7„) and depend on integer j. For all sets, pairs of parameters were again 
carefully chosen so that the actions of configurations with different j in the set were not 
much different. 























For APW configurations one expects: 


- log 


1 _ M 

Ni 


tot 


= RrP + const. = \(ot n + 7 n j) x oj(n) + const. 


(42) 


For a particular wave number n, one can plot [— log(iV^/iV tot )] vs. \(a n + 7 n j), and 


determine the slope uj(n) from a fit of the form (42). The expected linear dependence was 


seen in all our data sets at all couplings, wave numbers, and parameter choices; an example 
is displayed in figure [ll| 

Our main goal is to compare computed relative weights of test configurations with predictions 
of the DR ([ 5 ]) and GO (23) wave functionals. As already mentioned, NAC configurations are 
not suitable for that purpose, they only served us as test bed of our computer code. However, 
for APW the DR prediction for the dependence of the extracted function oj(n) on the wave 
number n: 

u(n) ~ fj, k 2 (n) ... dim. reduction (43) 


differs from our form: 


w(n) 


k 2 (n) 


yjk 2 (ri) + m 2 


GO proposal, 


where the wave momentum k above fulfils 


k 2 (n) = 2j2 1 - 


cos 


27rrij 

~L~ 


(44) 


(45) 


We therefore fitted our data for u{n) at each coupling /3 by the following two-parameter forms: 


uj{n) = 


a + bk 2 (n) 
ck 2 (n) 

yjk 2 (n) + m 2 


dim. reduction, 
GO proposal. 


(46) 


The result at (3 = 2.5 (30 4 lattice) is displayed in figure [l2| (green dashed line for DR, blue dotted 
one for our proposal). Both forms fit the data quite well at low plane-wave momenta, none of 
them is satisfactory for larger momenta. The situation at other gauge couplings is similar. 

Considerable improvement at all couplings is achieved by adding another parameter d to our 
form: 

ck 2 (n) 

yjk 2 (n ) + m 2 

see the red solid line in figure [T2| This would correspond, in the continuum limit, to adding 
a term 

d ( , „~ P2 .~ A ° n ) 7 (48) 


oj(n) = 


[1 + dk(n )], 


(47) 


{-V 2 - A 0 ) + m 2 


to the adjoint kernel ^ that appears in the VWF (23). 


For constant configurations with small amplitude the DR and GO forms coincide. For 
consistency between our data for NAC and APW configurations, the value of /xnac determined 
from sets of NAC configurations should agree with the appropriate combination (2cAPw/ r ^APw) 
of parameters obtained for abelian plane waves. Our results clearly pass this check quite success¬ 
fully, both values converge to each other at large /?, as exemplified in figure [13J 














P=2.5, 30 4 lattice 



Figure 12. oj(n) vs. /c(n) with fits of the 
forms (46) and (|47l) for APW configurations. 


Figure 13. The combination (2 c/m)f(/3) of 
the best fit (47) to APW data. Also displayed 
are values of /i/(/3) extracted from NAC data. 
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Figure 14. Parameter c (red) and rescaled Figure 15. Rescaled parameter df(/3 ) (green 
mass m/f(/3) (blue) of the best fit (47) vs. /3. points) of the best fit (47) vs. (5. 


Finally, the parameters of the best fit ©, c, m and d, if they correspond to some physical 
quantities in the continuum limit, should scale correctly when multiplied by the appropriate 
power of the asymptotic-freedom function /(/3) in (39). While the scaling of the ratio (2 c/m) 
multiplied by /(/?) has already been seen almost perfect in figure 13, it is not so convincing 
for individual parameters c and m/f(/3), though their growth in the range of f3 = 2.2 -4- 2.5 
is not large (see figure 14) and one may optimistically hope that it will level off at still higher 
couplings. On the other hand, the parameter d multiplied by f(/3) falls down quite rapidly in the 
same region (figure 15). The data thus suggest that the physical value of d might vanish in the 
continuum limit. This signals a possibility that the form (23) of the VWF might be recovered 
in the continuum limit. 


4. Open end of the romance; roads to take 

Let me summarize what has been achieved until now, both positive (^) and negative (^) 
points: 












































We have proposed an approximate form 
of the SU(2) YM vacuum wave functional 
that looks good in D = 2 + 1, somewhat 
worse in (3 + 1) dimensions. 

The method of Greensite and Iwasaki 
allows one to compute numerically (on 
a lattice) relative probabilities of various 
gauge-field configurations in the YM 
vacuum. 

For non-abelian constant configurations 
and for long-wavelength abelian plane 
waves the measured probabilities are con¬ 
sistent with the dimensional-reduction 
form, and the coefficients fi for these sets 
agree. 

Numerical data for test configurations 
are nicely described by a natural modifi¬ 
cation of our proposal, and the correction 
term seems to vanish in the continuum 
limit. 


In (3 + 1) dimensions, a method of gener¬ 
ating configurations distributed accord¬ 
ing to the proposed vacuum wave func¬ 
tional is not (yet?) available. 

This relative-weight method is only ap¬ 
plicable to compute weights of configura¬ 
tions rather close in configuration space. 


Neither the dimensional-reduction form, 
nor our proposal for the vacuum wave 
functional describe data satisfactorily for 
larger plane-wave momenta. 


No configurations tested so far, neither 
non-abelian constant nor abelian plane 
waves, can be considered typical. They 
are not in any sense true representatives 
of fields inhabiting the YM vacuum. 


The investigation reported in this talk could be extended in various ways: 

- One should compute, by the relative-weight method, weights of more realistic, “typical” 
configurations. One can e.g. generate ensembles of configurations of the full YM theory, take link 
matrices from a random time slice, make their Fourier decomposition, switch on/off or modify 
individual momentum modes, and compare the obtained momentum dependence of relative 
weights of configurations with that following from our, or some other, Ansatz for the VWF. 

- It is most desirable to find a way of generating field configurations distributed according 
to (the square of) the VWF. What was possible in D = 2 + 1 dimensions, is made much more 
complicated by the Bianchi constraint in (3 + 1). 

- If one achieved progress in the above points and gained more evidence for the proposed 
form of the vacuum wave functional in (3 + 1) dimensions, a number of other questions would 
call for an answer, e.g.: What are the dominant field configurations in the YM vacuum? Where 
do colour screening and IV-ality dependence arise from? ... 

One can just hope, on the crooked road towards understanding the QCD vacuum wave 
functional, not to run into a “Dead End” sign. 
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Appendix A. Bird’s eye view of the YM theory on the lattic^j 

In the compact lattice formulation, gauge fields A^(x) are represented by link matrices U^{x) 
(see figure [Al]): 

A m (x) = A“(x)T a (x) —» Uft(x) = exp[igaA M (x)]. (A.l) 

The field strengths F /ti ,(x) are related to products t/ MI ,(x) of link matrices along a plaquette: 

F^(x) = F^(x)T a (x) —► U^{x) = U^x)U v {x + jl)Ul{x + u)Ul{x). (A.2) 



Figure Al. Building blocks 
of pure Yang-Mills theory on 
a lattice. 


The simplest form that replaces the euclidean SU(A) YM action 

S = t / d 4 x Tr [F IJiV (x)F ^(x)} 


(A.3) 


is the Wilson action 


S w = a 4 pY, 1 

p L 


1 

N 


Re Tr U P 


with (3 = 2 N/g 2 . 


(A.4) 


Vacuum expectation values 


<0|Q|0> = J{dU}%{U}Q* 0 [U] (A.5) 

are given by path integrals 

(Q) = tf /Q[*7] exp(-5w[^])- (A.6) 

In the numerical Monte Carlo simulation one computes in fact: 

I -^conf 

<«> * ^ £ <2[{C.}], (A.7) 

an average over a (large) number V con f of gauge-field configurations {Cj} distributed according 
to the probability distribution ~ exp(—5^1^]) ~ |\D , o[L]| 2 - 


Appendix B. Recursion method for simulation of the vacuum wave functional 

Let us define a probability distribution for gauge fields A in a background of a second, 
independent configuration A': 


V[A-R[A']} = J\f [A 7 ] exp 


d 2 x cPy B a (x- A) J^A'] B b (y ; A) 


(B.l) 


7 For more details see e.g. the textbook (56] . 
















where A and A' are fixed to a variant of axial gauge (Appendix B of |23j). If we assume that 
the variance of H is small among thermalized configurations, we can approximate: 


P[ A] = V[A-,&[A}} « V[A- (£)] « j dA ' V[A; &[A']] P[A']. 


Then the probability distribution can be generated by solving (B.2) iteratively: 


P«[A] = P[A;£[0]], 

P (fc+1) [A] = J dA'T[A-,&[A'}] pW[A']. 


(B.2) 


(B.3) 

(B.4) 


A block diagram of a practical implementation of this recursion procedure is shown below: 



In the case of the kernel appearing in the vacuum wave functional (12), the procedure converges 
quite rapidly, typically in 0(10) above cycles. The hypothesis about small changes of R among 
equilibrated configurations is confirmed a posteriori by the absence of large fluctuations of the 
spectrum of for individual recursion lattices. 
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